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We give results for B, B s , B c and bottomonium spectroscopy using NRQCD heavy quarks and 
HISQ valence and sea quarks. Five MILC ensembles of gluon configurations with three values of 
the lattice spacing and miight/'Wstnmge values down to 0.1 are used that include 2+1+1 flavours of 
sea quark. Systematic errors in the NRQCD action are improved through the radiative correction 
of the coefficients of terms at v 4 . Improved results for S-wave and P-wave bottomonium states are 
discussed as well as a prediction for the full D-wave spectrum. Preliminary results for the ground 
state B meson masses are also presented. 
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1. Introduction 



The spectroscopy of mesons containing b quarks has been an important testing ground for 
lattice QCD [1] and recent improvements in a number of areas mean it is necessary to revisit 
the bottomonium and B meson spectra. The study of the bottomonium spectrum on the lattice is 
simpler than for lighter hadrons in several important ways. There are a number of states below 
the threshold for decay into two B mesons meaning that for low lying states there is no systematic 
error from not including these multiparticle diagrams on the lattice. The splittings between excited 
states in the spectrum show little dependence on the valence b quark mass, which reduces the error 
due to mistuning, or the light sea quark masses, which reduces errors from chiral extrapolation. 

The use of effective field theories such as NRQCD for b quarks avoids the large discretisation 
eiTors that are typical of other lattice quark formulations on coarse lattices. NRQCD is also com- 
putationally very cheap since propagators can be calculated from a time evolution rather than an 
inversion of the Dirac matrix. In recent years the HPQCD collaboration have calculated order a s 
radiative corrections to the Wilson coefficients in the NRQCD action which were the main source 
of systematic error. We present a new study of the bottomonium spectrum, along with preliminary 
results for the lowest lying B meson states, using these coefficients and a number of other improve- 
ments over the previous HPQCD calculation. This serves as a check of the improvements made 
and allows the remaining parameters, such as the mass and lattice spacing, to be fixed. We are 
also able to make a prediction for the D-wave states. Successfully reproducing the known part of 
the spectrum gives us confidence in more phenomenologically interesting quantities, such as decay 
constants and mixing matrix elements, which will follow on from this calculation. 

Further improvements have been made in the gluon and sea quark sector. We use the recent 
MILC collaboration ensembles with 2+1+1 flavours of HISQ sea quarks at three lattice spacings 
(from 0.1 5fm to 0.09fm) and two light quark masses (mi/m s = 0.1,0.2) [2]. These have more 
chiral light quarks, include charm for the first time and have a larger spatial volume than previous 
ASQTAD ensembles. The coefficients of the gluon action have also been perturbatively improved 
[3], including the effect of the HISQ sea quarks. 

Full details of the S and P-wave results presented in these proceedings can be found in [4]. 



2. Calculation details 

Heavy quark action: NRQCD is an expansion of QCD in powers of the heavy quark velocity v. 
The form of the NRQCD Hamiltonian is: 

A (2) 
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This is the standard action including terms of order v 4 and some other terms designed to reduce 
discretisation errors, see ref. [5]. am,}, is the b quark mass in lattice units, n is the stability parameter 
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used in the evolution equation, E,B are the improved chromo-electric and chromo-magnetic field 
strengths and the AW are lattice derivatives. These are all defined in more detail in [1]. With the 
coefficients c, = 1, this action is equivalent to QCD at tree level. The c, can be expanded in powers 
of a s , as Ci = 1 + cj ■ cc s + ^{ccj) and most of the have now been calculated perturbatively 
by HPQCD by matching to QCD at one loop. Radiative corrections are kept small by tadpole 
improving the gauge fields in the action with the Landau link mol- Tuning of some of the coefficients 
has also been performed nonperturbatively as an independent check and to estimate higher order 
contributions. Radiative corrections are particularly important for quantities such as the hyperfine 
splitting, which is proportional to c\, and were the dominant source of systematic error in previous 
calculations. 

Tuning of parameters : The lattice spacings were fixed to a precision of better than 1 % using the 
Upsilon 2S-1S splitting. For a consistency check, the static quark potential parameter r\ was also 
calculated using the decay constant of the tj v meson [6] and MILC values for r\ja. Both methods 
give continuum values for r\ that agree. 

The b quark mass was tuned by computing the spin averaged kinetic mass = I (3Mj + M Tlb ) 
and comparing to the experimental value appropriately adjusted for missing electromagnetic and 
annihilation effects. In ref [4], we present a detailed study of systematic errors in the tuning of 
amb and provide accurate values for each ensemble. The light, strange and charm quarks in the 
heavy-light calculations were included with the highly improved staggered quark (HISQ) action 
[7] and the strange and charm quark masses were tuned to the masses of the T] v and tj c . Stochastic 
noise sources were used for the S and P-wave propagators and all B mesons. Data was extracted 
from the correlators using a multi-exponential matrix Bayesian fit. Full details of all parameters 
used are given in [4] . 



Here we present a subset of the results obtained in [4] starting with an overview of the bottomo- 
nium spectrum obtained in figure 1. Not all points on the plot are pre/post-dictions, the Upsilon 
2S — IS splitting was used to fix the lattice spacing and the spin average of the IS states was used 
to tune amb- 

Splitting ratios: The quantities which can be calculated most accurately are ratios of energy 
splittings. Some systematic errors cancel in the ratio leaving an estimated systematic on each 
ensemble of < 1%. Results are plotted in figure 2, and we extract physical values of 



which are in good agreement with experiment (from [8]) and have overall errors of 2.4% and 1.4% 
respectively. 



3. Results 



T(35-15) 
T(25-15) 
l l Pi -IS 
T(2S-15) 



0.820(12), Expmt = 0.8088(23) 



1.625(39), Expmt = 1.5896(12) 



(3.1) 
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Figure 1: An overview of the spectrum obtained 
in our calculation. The T(1S,2S) are used for fix- 
ing the lattice spacing and tuning the b quark mass. 
Other points are pre/post-dictions. 
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Figure 2: Plot of the ratios of energy splittings for 
each ensemble along with the physical value ob- 
tained from the fit. 



P-wave splittings and nonperturbative tuning: The P-wave spin splittings are shown in figure 
5. We tune £3,04 (see [1]) so that these splittings agree with the experimental values (crosses). The 
plot also shows the tree level coefficients which give a 3 Po — 3 P splitting that is too small. This 
nonperturbatively tuned C4 agrees with the perturbative calculation within errors. The fine lattice 
P-wave masses are shown on figure 1 for the perturbative value of C4. 

Hyperfine splittings: Another important test of any computation of the spectrum is the hyper- 
fine splitting M(T) —M(r\b). Previous tree level lattice NRQCD calculations suffered from large 
systematic errors of around 25%. The dominant contributions to this error come from radiative cor- 
rections, since the splitting is proportional to c|, and &(y k ) terms in the action. Figure 3 shows our 
new results for the hyperfine splitting. The plot shows the result from both perturbatively and non- 
perturbatively (using the experimental P-wave spin splittings) calculated C4 values and the physical 
result extracted from this data. The values have been adjusted for mass mistuning and the effect 
of missing 4-quark operators has been calculated perturbatively. Correlated systematic errors for 
missing higher order terms in a s were included in the fit. Finally, a 10% systematic error has been 
applied for missing v 6 terms. 

The ratio of the 25 to 15 hyperfines is shown in figure 4. This ratio should be independent of 
C4, a fact which is confirmed by the data, so the error is dominated by statistics. Our results for the 
15 and 25 hyperfine splittings are: 

M(T(15)) -M{ri h (lS)) = 70(9)MeV, M(T(25)) -M(rj b (2S)) = 35(3)MeV 

D-wave splittings: We are also able to expand upon previous NRQCD spectrum calculations by 
giving a prediction of the D-wave spin splittings. Several different lattice operators were used for 
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Figure 3: Results for the hyperfine splitting. Pertur- 
bative and nonperturbative C4 values are shown with 
errors from statistics, lattice spacing, mass retun- 
ing and 4-quark operators shown. Correlated higher 
order systematic errors are not shown on the data. 
Tree level C4 = 1 points are also plotted but are not 
included in the fit 



Figure 4: Results for the ratio of 2S to IS hyperfine 
splittings from a 5 x 5 matrix fit. Pink open squares 
show the result for a tree level C4 and blue crosses 
denote results using the nonperturbatively tuned val- 
ues. The corresponding ratio for charmonium is in- 
cluded for comparison. 



each physical state and the result for the coarse ensemble is shown in figure 6. The splittings have 
no noticeable lattice spacing or light quark mass dependence. 

The ratio mi,/m s : Using the tuned values of ami, an d am s we can calculate the ratio of the two 
quark masses. These are converted to the MS scheme at scale /I via the pole mass. The NRQCD and 
HISQ mass renormalisations are known to one loop and the dominant error comes from missing 
&(ccf) terms. The result obtained is 54.7(2.5) and is shown in figure 7. Also included is a previous, 
independent HPQCD value (53.4(9)) using HISQ for both the b and s quarks [9, 10]. 

B meson hyperfine splittings (preliminary): Here we discuss preliminary results for the low 
lying B meson states. The B meson masses are extracted from NRQCD-HISQ correlators with 
random noise sources using local and exponential source smearings with two different radii. In 
figure 8 we provide further evidence to support the previous HPQCD prediction in [11] that ratios 
of B meson hyperfine splittings are equal to one - as above, the ratio should be independent of C4. 
This implies that the heavy-light hyperfine splitting is independent of the light valence quark mass. 

4. Discussion 

The calculation of the bottomonium spectrum presented here shows a significant improvement 
in both systematic and statistical uncertainty over the previous HPQCD results with evidence for 
similar levels of improvement to come in the heavy-light sector. This is part of a broader heavy 
quark physics program in progress by the HPQCD collaboration and gives us the opportunity to 



5 



Bottomonium spectroscopy from lattice NRQCD 



R. J. Dowdall 




4± 



set 3 c,=l 



Figure 5: Results for the P-wave spin splittings rel- 
ative to the spin averaged i P state. Set 3 and 4 label 
the two coarse ensembles (~ 0.12fm) in our calcu- 
lation with mi/m s = 0.2,0.1 respectively. The two 
values for the 3 f2 state denote the E and T2 lattice 
irreps (left and right respectively.) 
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Figure 6: Plot of the D-wave bottomonium states 
relative to the spin averaged 3 D state on a coarse 
ensemble. Errors are dominated by statistics. 
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Figure 7: Plot of mi,/m s in the MS scheme at scale 
11. The errors on the points include statistical, lat- 
tice spacing and NRQCD systematics. The previous 
HISQ result is shown with a black circle. 
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Figure 8: A plot of the ratio (Mg* — M Bs ) / (M B * - 
Mg.) (where q = d,c) on two coarse (set 3,4, a ~ 
0.12fm) and a very coarse (set 1, a ~ 0.15fm) en- 
semble. The previous HPQCD value of 1 .00(23) is 
included for comparison. 
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accurately tune all necessary parameters and perform checks of the perturbative improvements. 
The main aims of the program are improved determinations of important quantities for B physics 
phenomenology using the NRQCD and HISQ actions. Our results demonstrate that our errors are 
well understood and that subsequent computations will be reliable. 
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